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^ : ABSTRACT 

■ We report observation results of the prompt X- and 7-ray emission from 

l/^ ' GRB 000528. This event was detected with the Gamma- Ray Burst Monitor 

and one of the two Wide Field Cameras aboard the BeppoSAX satellite. The 

GRB was promptly followed on with the BeppoSAX Narrow Field Instruments 

and with ground optical and radio telescopes. The X-ray afterglow source was 

t~^ ■ identified, but no optical or radio counterpart found. We report here results 

^ ■ from the prompt and afterglow emission analysis. The main feature derived from 

O 

O ■ density with evidence of its decrease with time. We model this behavior in terms 

rS^ • of a time- dependent photoionization of the local circum-burst medium, finding 

I ■ that a compact and dense environment is required by the data. We also find a 

^ ■ fading of the late part of the 2-10 keV prompt emission which is consistent with 

^ ■ afterglow emission. We discuss this result in the light of the external shock model 

cz ■ 

scenario. 
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1. Introduction 

The prompt emission of long duration (>2 s) Gamma Ray Bursts (GRBs) in the broad 
2-700 keV energy band, which was accessible to the BeppoSAX wide field instruments (Wide 
Field Cameras, Jager et al. 1997 , and Gamma Ray Burst Monitor, Frontera et al. 1997), 
has been demonstrated to provide information not only on the X-ray emission mechanisms 
of GRBs but also on their environments (e.g., Frontera et al. 2000, Frontera 2004a). The 
properties of these environments are of key importance to unveil the nature of the GRB 
progenitors, which is a still open issue. Many observational results point to the collapse 
of massive stars, in particular to the explosion of type Ic supernovae, as the most likely 
progenitors of long GRBs. However the simultaneity of the GRB event and supernova 
explosion (hypernova model, e.g., Paczynski 1998), in spite of the convincing results obtained 
from GRB0303289 (Stanek et al. 2003; Hijorth et al. 2003), cannot be considered a solved 
problem for all GRBs. A two-step process, in which first the massive star gives rise to a 
supernova explosion with the formation of a neutron star, and thus, after some time, to 
a GRB from the delayed collapse of the NS to a BH (supranova model, Vietri & Stella 
1998) or from a phase transition to quark matter (Berezhiani et al. 2003) cannot be robustly 
excluded. At least in some cases (e.g., GRB990705, Amati et al. 2000 ; GRB011211, Frontera 
et al. 2004b, GRB991216, Piro et al. 2000), in which X-ray transient absorption features or 
emission lines have been observed and found consistent with an Iron-enriched environment 
typical of a previous supernova explosion, the two-step process is favored (Lazzati et al. 
2001a, 2001b, 2002). 

The properties of the GRB environment can also be derived from the study of the 
continuum spectrum of the prompt X-ray emission, which is expected to show a low-energy 
cutoff in case of a high column density along the line of sight. If the absorbing material is 
located between a fraction of a parsec to several parsecs from the GRB site, the cutoff is 
expected to fade as the GRB evolves. This is due to the progressive photoionization of the 
gas by the GRB photons (e.g., Bottcher et al. 1999; Lazzati & Perna 2002). This is an effect 
analogous to the one originally proposed by Perna & Loeb (1998) in the optical domain with 
resonant lines. Actually, such cutoffs, with also evidence of their decreasing behaviour with 
time from the GRB onset, have already been detected (GRB980506, Connors & Hueter 1998; 
GRB980329, Frontera et al. 2000). Lazzati & Perna (2002) interpreted them as evidence 
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for presence of overdense regions in molecular clouds with properties similar to those of star 
formation globules. 

Now we report the discovery of another case of absorption. It has been found in the X- 
ray spectrum of the prompt emission of GRB000528, in the context of a systematic investiga- 
tion (Frontera et al. 2004c) performed to study the spectral evolution of all GRBs jointly de- 
tected with WFCs and GRBM. GRB000528, detected with the GRBM (Guidorzi et al. 2000) 
and localised with the BeppoSAX WFC (Gandolfi 2000), was promptly followed- up with the 
BeppoSAX Narrow Field Instruments (NFI), with the discovery in the WFC error box of a 
likely X-ray afterglow candidate of the burst (Kuulkers et al. 2000). However neither an op- 
tical or radio counterpart was found (lowest upper limit to the R magnitude of 23.3, Palazzi 
et al. 2000; 3.5(7 upper limit of 0.14 mJy at 8.46 GHz, Berger & Frail 2000). 

In this paper we report our findings and their interpretation. We also discuss the follow- 
up X-ray observations performed with the BeppoSAX NFIs and their implications. 



2. BeppoSAX observations 

GRB000528 was detected with the BeppoSAX WFC No. 2 and GRBM (Gandolfi 2000) 
on 2000 May 28 at 08:46:21 UT. The GRB position was determined with an error radius of 
2' (99% confidence level) and was centered at aaooo = 10''45"^06.3^ and (52ooo = -33°58'59" 
(in 't Zand et al. 2000). The burst was then foUowed-up with the NFIs two times, the first 
time (TOOl) from May 28.71 UT to May 29.57 UT (from 8.25 to 28.90 hrs after the main 
event), and the second time (T002) from May 31.65 UT to June 1.49 UT (from 3.284 to 4.124 
days from the main event). Two previously unknown sources, the first one (SI) designated 
ISAX J1045.5-3358, with coordinates ^2000 = 10M5'"28% and ^2000 = -33°58'14", the 
second one (S2) designated ISAX J1045. 1-3358, with coordinates aaooo = 10^45'^08^ and 
(^2000 = — 33°59'26" were detected with MECS and both found to have disappeared in T002 
(Kuulkers et al. 2000). However only the S2 position (uncertainty radius of 1.5') was found 
consistent with the WFC error box (see also Hurley et al. 2000). 

We have reanalysed the MECS and LEGS data. The exposure times of the two instru- 
ments to the GRB direction are 8084 s (LEGS) and 26589 s (MECS) in TOOl, and 6905 s 
(LEGS) and 31210 s (MECS) in T002. The images of the sky field derived with MECS 2+3 
during the two TOOs are shown in Fig. 1, with superposed the WFC error box. The two 
sources reported by Kuulkers et al. are clearly visible in the left image. The probability that 
the excesses found are due to background fluctuations is ~ 7 x 10^^. The source SI, which is 
out of the WFC error box, is coincident with a radio source, which was already present in the 
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NRAO VLA Sky Survey (NVSS) at the time of the GRB occurrence. It was again detected 
by Berger & Frail (2000) during the search of the radio counterpart of GRB000528. This 
source is also detected in the second observation (chance probability of ~ 5 x 10~^), with 
a flux variation by a factor 1.7 from the first to the second TOO. However the source S2 is 
detected only during the first TOO, it is consistent with the WFC position, and undergows 
a flux decrease larger than a factor 3 from the first to the second TOO (see Section 3.2). As 
already concluded by Kuulkers et al. (2000), we assume this object as the likely afterglow 
source of GRB000528. 



3. Data analysis and results 

The spectral analysis of both the prompt and afterglow emission was performed with 
the XSPEC software package (Arnaud 1996), version 11.2. All the reported parameter uncer- 
tainties are given at 90% confidence level. 



3.1. Prompt emission 

Data available from GRBM include 1 s ratemeters in the two energy channels (40- 
700 keV and >100 keV), 128 s count spectra (40-700 keV, 225 channels) and high time 
resolution data (up to 0.5 ms) in the 40-700 keV energy band. WFCs were operated in 
normal mode with 31 channels in 2-28 keV and 0.5 ms time resolution. The burst direction 
was offset by 8.0° with respect to the WFC axis. With this offset, the effective area exposed 
to the GRB was ~ 340 cm^ in the 40-700 keV band and 70 cm^ in the 2-28 keV energy band. 
The background in the GRBM energy band was estimated by linear interpolation using the 
250 s count rate data before and after the burst. The WFC spectra were extracted through 
the Iterative Removal Of Sources procedure (IROS ^, e.g. Jager et al. 1997 ) which implicitly 
subtracts the contribution of the background and of other point sources in the field of view. 

Figure 2 shows the light curve of GRB000528 in two energy bands. In X-rays (top 
panel of Fig. 2) it shows a seemingly long rise with superposed various small pulses until 
it achieves the highest peak at a level of 300 cts s~^ after ~ 70 s from the onset. Then 
the X-ray flux starts to smoothly decay and ends after ~ 120 s from the flux peak. In 
gamma-rays the light curve (bottom panel of Fig. 2) is totally different. In correspondence 
of the first of the secondary pulses of the X-ray light curve we find the most prominent 
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gamma-ray pulse, where the GRB achieves the maximum count rate in the 40-700 keV band 
of 1700 counts s~^. Instead, in correspondence of the peak flux in X-rays we only find a 
secondary pulse of the gamma-ray light curve. Smaller and shorter duration pulses are also 
visible in the gamma-ray light curve. However the smooth decay of the X-ray light curve is 
not observed in gamma-rays and the entire duration of the 40-700 keV burst is only ~ 80 s 
versus the ~ 160 s duration in X-rays. 

For the study of the GRB spectral shape and its evolution with time, we extracted the 
WFC and GRBM data in the same time intervals and we simultaneously fitted the WFC 
and GRBM data. For GRBM we used both types of data available: the 1 s counters and 
the 128 s spectra. Using the counters we had the possibility to subdivide the GRB light 
curve in 6 adjacent time intervals A, B, C, D, E, F of 18, 11, 10, 20, 20 and 90 s duration, 
respectively (see Fig. 2) and to derive 2-700 keV count spectrum for each of these intervals. 
The obtained GRB count spectra are shown in Fig. 3. Given that the first three spectra (A, 
B, and C) have a low statistical quality in the 2-28 keV band, for the fitting analysis they 
were summed together (interval ABC). Using 128 s GRBM transmitted spectra, we obtained 
two 2-700 keV spectra in the two time intervals 1 and 2, whose boundary time is marked in 
Fig. 2. The interval 1 covers the ABC interval above, while the interval 2 covers the intervals 
D, E and F. 

Each of the four spectra (ABC, D, E, F) was fit with various continuum models: a 
power-law (pl), a power-law with an exponential cutoff (cpl), a broken power-law (bknpl), 
a smoothly broken power-law (or Band law, BL, Band et al. 1993). All these models were 
assumed to be photoelectrically absorbed by an equivalent-hydrogen column density Nn 
to be determined (the WARS model in xspec). In the fits, the normalisation factor of the 
GRBM data with respect to WFC was left free to vary, checking that its best fit value was in 
the range (0.8 to 1.3) in which this parameter was found in the fiight calibrations with either 
celestial sources (e.g.. Crab Nebula) or GRBs observed with both the WFC-GRBM and 
the BATSE experiment aboard the Compton Gamma Ray Observatory (e.g., Fishman et al. 
1994). The power-law model gave a completely unacceptable fit of all count spectra either 
assuming the Galactic column density along the GRB line of sight (A^^ = 6.1 x 10^" cm~^) or 
leaving the N-^ free to vary in the fit. The fit with the Band law gave in general unconstrained 
parameters. The other two models (cpl), and bknpl) resulted the most suitable to describe 
the data. The fit results with these models are reported in Table 1, along with the x^/dof 
value either in the case Nn is frozen to the Galactic value (A^h ) o^ ^^ ^^e case it is left free 
to vary in the fit. 

Apart from the typical hard-to-soft evolution of the spectrum with time, the most re- 
markable feature found is the value of the column density Ah, which results to be higher 
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than the Galactic value in the first three time intervals ABC, D and E (see chance proba- 
bilities obtained with F-test in Table 1) and it is consistent with 0, and thus with N^, in 
the interval F {2a upper limit reported in Table 1). As can be seen from the F test values, 
these results are independent of the spectral model chosen (bknpl or cpl). We notice that, 
altough with lower significance, these results are confirmed with fits using WFC data only. 
Also the analysis of the two GRBM 128 s spectra are consistent with these results. Both 
spectra were well fit with either a bknpl or cpl. In both cases a column density A^h = A^h 
does not give the best description of the data (see Table 1); a best fit is obtained with A'h 
values significantly higher than the Galactic column density. 

Even if we cannot completely exclude a constantly high A^h, the data are in favour 
of a decreasing column density. Indeed, summing together the data in the intervals ABC 
and D, we estimate an A^h = l-O^o's x 10^^ cm~^, which is not consistent with the highest 
2o" upper limit of the A'h in the interval F (4.4 x 10^^ cm~^) and it is only marginally 
consistent with its Scr upper limit. In addition, we find that the Ah centroid is continuously 
decreasing with time (see Fig. 4) and the fit of the A'h values with a constant (weighted 
mean value 5.8 ± 1.1 x 10^^ cm~^) gives a x^/dof value of 7.8/3, versus a value of 0.65/2 in 
the case of an exponential law (A'H(t) = Aoexp(— t/to)). Using the F-test (Bevington 1969), 
the probability that a constant A'h value is the real parent distribution is 6%. Finally the 
comparison of the light curves in the 2-10 keV and 10-28 keV energy bands (Fig. 4) are also 
suggestive of a variable A^h- Assuming the above exponential law, the best fit parameters of 
the exponential law are Aq = (2-4 ±1.2) x 10^^ cm~^ and to = 51 ± 23 s. 

No evidence of an absorption edge associated with the variable A^h is found in the 
spectra. The upper limit to the edge optical depth r depends on the energy assumed and on 
the time interval considered. With the continuum parameters fixed to their best fit values 
obtained with the CPL, the 2a r upper limit achieves the lowest value (0.2) at 5 keV in 
interval D and the highest value (1.8) at 2 keV in the interval ABC. The other upper limits 
are in between. 

On the basis of the best fit results, peak fiuxes and fiuences associated with the event 
were derived. The 2-10 keV and 40-700 keV peak fiuxes in 1 s time bin are -F2-10 = 2.1 x 
10"^ erg cm~^ s~^and -F40-700 = 1-4 x 10~^ erg cm~^ s~^, while the fiuences are given by 
5'2-io = 9.1 X 10^'^ erg cm^^, and 5*40-. 700 = 1-37 x 10^^ erg cm^^, with a softness ratio SR = 
5'2-io/'S'4o-7oo = 0.066, which is a value typical of a classic GRB (e.g., Frontera et al. 2000). 
Finally the total 2-700 keV fiuence of the GRB is given by 6*2-700 = 1-9 x 10~^ erg cm~^. 
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3.2. Afterglow emission 

Due to the weakness of the afterglow source (see Fig. 1), no statistically significant 
spectral information could be derived. Assuming the average power-law photon index (1.95) 
derived from the spectra of a large sample of GRB afterglows (Frontera 2003a), from the 
MEGS data we estimated the afterglow source fiux in the two observations of GRB000528, 
using an image extraction radius of 2' and adopting the MEGS standard background level. 
The resulting 2-10 keV mean fiux is (1.6 ±0.4) x 10^^'^ erg cm^^ s^^(lo" uncertainty) during 
TOOl and less than 5.2 x 10^^^ erg cm^^ s^^(2o' upper limit) for T002. These fiuxes can 
only give un upper limit to the fading rate of the afterglow source. Assuming a power law for 
the fiux time behaviour {F(t) oc t^), with time origin at the GRB onset, the 2a upper limit to 
6 is —0.74, which is not very constraining for the afterglow decay rate. Back extrapolated to 
the time of the main event, this decay implies a fiux much lower than that measured with the 
WFG. A tighter constraint on 6 can be obtained with a reversed argument, imposing that - 
as observed in many BeppoSAX GRBs (Frontera et al. 2000; Frontera 2003a) - the late part 
of the prompt emission is correlated with the late afterglow emission, pointing to the former 
as early afterglow. In Fig. 5 we show the results. Assuming the GRB onset as origin of the 
afterglow onset time, we find that the smooth GRB decay starting at the beginning of the 
time interval F in Fig. 2 is consistent with a power-law with index 6 = —3.7 ±0.3. As shown 
in Fig. 5 (left panel), its extrapolation to the time of the BeppoSAX NFI measurements gives 
a fiux much lower than that measured. A different result is obtained if we fix as origin of the 
time coordinate for the afterglow the beginning of the time interval F. In this case, the 2-10 
keV smooth decaying fiux is described by a broken power-law (a simple power-law gives 
is unacceptable, x^/dof = 88/31, see right-hand panel of Fig. 5), with best fit parameters 
6i = —0.12 ± 0.06, 62 = —1.33 ± 0.13 and break time threak = 16l2;5 s- ^^ notice that the 
second power-law index is that typical of X-ray afterglows (Frontera 2003a) and that, when 
this law is extrapolated to the epochs of the TOOl and T002 measurements, it is found 
consistent with the afterglow fiuxes measured at these times. All that is suggestive of the 
fact that the afterglow onset time is delayed by ~ 80 s with respect to the GRB onset time 
and that the final fading of the early afterglow of GRB000528 occurs about 16 s after a 
slower fading. 



4. Discussion 

Two main results have been obtained from the present analysis of the GRB000528 
BeppoSAX data: the evidence of an early afterglow with onset time delayed with respect to 
the GRB onset, and the discovery of an absorption column density significantly higher than 



the Galactic value during most of the burst prompt emission with evidence of a decreasing 
behavior with time. We discuss both. 



4.1. The early afterglow 

The finding that, assuming as time origin the beginning of interval F in Fig. 2, the 
extrapolation of fading law of the late prompt emission is consistent with the flux of the 
late afterglow, is highly suggestive of the expectations of the internal-external shocks model, 
during the fast cooling phase (Sari & Piran 1995; Sari & Piran 1999). The above derived 
early afterglow decay law and the photon spectrum of the last interval F (see Table 1) confirm 
this scenario. Indeed, assuming a synchrotron emission, according to Sari & Piran (1999), 
the photon spectrum due to the fast cooling electrons is described by a set of four power- 
laws, one N{E) (x v ioi v < v < Vsa-, where Vga is the photon self-absorption frequency, 
the second N{E) oc z/~^/^ for u^a < t^ < t'o where u^ is the cooling frequency, the third 
N(E) oc z/^'^/^ for z/c < z/ < Um, where z/^ is the characteristic synchrotron frequency, and 
the last N{E) oc z/~^/^~^ for z/ > z/^, where p is the index of the electron energy power- 
law distribution. Comparing these expectations with our spectral results in the interval F 
(broken power-law with low-energy photon index a = 1.27^o;33, high-energy photon index 
P = 4.6_3.i, Ehreak = 51^29); we cau See that a is consistent with the expected power law 
index (3/2) in the Uc < i' < i^m frequency band, and thus that z/^ corresponds to Ei,reak- 
Unfortunately the high-energy photon index f3 is affected by a large uncertainty with the 
90% upper boundary unconstrained, and thus it cannot be efficiently used to constrain p. 
However p can be constrained by the measured light curve of the early afterglow. Indeed, in 
the case of the fast cooling phase (Sari & Piran 1999), at X-ray frequencies, the energy flux 
F{t) is expected to increase as t^ for t < t^, where t^ is the time in which the relativistic flow 
changes from a constant Lorentz factor into a decelarating phase, while F{t) is expected 
to decay as t~^/^ for t^ < t < tm, where tm is the time in which the typical synchrotron 
frequency cross the observed frequency; finally F{t) is expected to decay as t^/^^^p/^ at 
later times (t > tm)- The light curve we observe is nicely similar to that expected: it has 
first a temporal index 6i = —0.12 ± 0.06, which is marginally consistent with that (-1/4) 
expected for t < tm, then has a temporal index 62 = —1.33 ± 0.13, which, compared with 
that expected (1/2 — 3p/4), allows to estimate p = 2.44 ± 0.17. This value of p is consistent 
with the estimated high-energy photon index of the bknpl. 

Due to the thick shell regime, we are not capable to directly observe the onset time of 
the afterglow and then to derive the initial Lorentz factor Fq from the decelaration time of 
the fireball t^, which depends on Fq (Sari & Piran 1999). An attempt to evaluate t^ was 
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done by fitting the GRB light curve in the interval F with a broken power-law as above, 
in which this time 6i is frozen to —1/4 and the other parameters, 62, tbreak and the onset 
time to of the afterglow with respect to the starting time of the interval F, were left free 
to vary. We found that the best fit is obtained with to = —3^5 s, S2 = —1.29to\l, ^^^ 
tbreak = 17^i s. The fading law obtained with these parameters is still nicely consistent with 
the late afterglow measured fiux. From this result, following Sari & Piran (1999), we can get 
an estimate of the deceleration time t^ (t^ = —to), and then an estimate of Fo if we know 
the GRB redshift. Placing GRB 000528 at z = 0.5 (see below for a discussion), we obtain 
Fo ~ 100n4 where n4 is the density of the external medium in units of 10'^ cm~'^. We 
stress however that the derived value of these parameters should be taken as rough estimates 
also because we adopted equations valid for an unperturbed interstellar medium, which are 
not strictly applicable to the thick shell case. In addition, it is worth considering that, 
lacking a decay slope only constrained from NFI observations, other possible solutions could 
be considered. In spite of all this uncertainty we consider that outlined a plausible scenario. 



4.2. The column density behavior 

The other major result of the present spectral analysis of the GRB000528 prompt emis- 
sion is the detection of an hydrogen-equivalent column density significantly higher than the 
Galactic value. Even if, as discussed, a constant value cannot be excluded, the likely tem- 
poral behavior of A^h is an exponential function with initial value A^h(O) ~ 2.4 x 10^^ cm~^ 
(assuming a redshift z = 0) and decay constant to ~ 51 s. It is the second time, after the 
case of GRB980329 (Frontera et al. 2000; Lazzati & Perna 2002), that A^h such time behav- 
ior, although evidences of variable absorption during the prompt emission of GRBs have 
already been reported (GRB990705, Amati et al. 2000, GRB010222, in 't Zand et al. 2001, 
GRB010214, Guidorzi et al. 2003). As discussed in Section 1, in the case of GRB980329, 
the A^H behavior was found consistent with a photoionization process, which is expected 
when a GRB occurs within a cloud of initially cold gas. The cloud mass density was found 
to be n ~ 4.5 x 10^ cm~^ for a composition typical of interstellar matter, and the radial 
distribution was consistent either with a uniform sphere of radius Rsphere = 0.13 pc, or a 
shell at distance Rsheii = 0.066 pc and width AR = O.lRgheii- A marginal but not conclusive 
preference was given to the shell geometry (Lazzati & Perna 2002). 

For GRB000528 we tested the same model with two important improvements in the 
simulation of the photoionization process. First, we have taken into account the true gamma- 
ray light curve of the GRB, while in the case of GRB980329 a constant fiux was assumed to 
be incident on the cloud. Second, we have taken into account the cosmological effect on both 
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the time dilation of the light curves and the measured column density: the fit was performed 
for three different values of z (0.1, 0.3, 0.5). Consequently, the A/'h values shown in Fig. 4 
were corrected to take into account the redshift effect on the column density according to 
the relation A^h(-2) = ^h(O) (1 + -z)^'^, which is valid in the redshift range < z < 4 
(Lazzati & Perna 2002). The results of the model fit to the data are shown in Fig. 6 and in 
Table 2, for the three values of z. From the x^ values shown in Table 2, also in this case 
the two cloud geometries (sphere or shell) cannot be significantly disentangled, even if it 
appears, also in this case, a slight preference for the shell geometry. Assuming a spherical 
geometry, the radius of the cloud increases with z, from 8.9 x 10^^ cm to 5.0 x 10^^ cm, while 
in the case of a shell, the shell distance increases with z from 5.6 x 10^^ cm to 1.8 x 10^^ cm 
(see Table 2). Also the initial column density is sensitive to the redshift assumed, ranging 
from from 6.3 x 10^^ cm~^ to 1.0 x 10^^ cm~^ in the first case, and from 3.5 x 10^^ cm~^ 
to 8.9 X 10^^ cm~^ in the second case^. However, independently of the redshift and cloud 
geometry, the best fit values are all consistent with overdense regions of molecular clouds 
where star formation takes place, as also found for GRB980329 (Lazzati & Perna 2002). 

The amount of mass required to explain the variable absorption ranges between few to 
few hundreds solar masses. This mass range is consistent with the masses ejected by Wolf- 
Rayet stars, the candidate progenitors of GRB explosions. If the mass were distributed as 
a wind, with an r~^ profile, most of the absorption would be at small radii, and the column 
density would be expected to vanish quickly. Lazzati & Perna (2002) simulated the evolution 
of absorption from a stellar wind, finding that even a Thompson thick wind would vanish 
in less than a second (see their Fig. 6). This does not mean that our observations provides 
evidence against a stellar wind. It means that this kind of observation is not sensitive to 
the presence of a wind. Nevertheless, the observation can be interpreted as a circumstantial 
evidence for the stellar origin of the progenitor. As recently discussed by Chevalier et al. 
(2004), the interaction of the wind with the surrounding ISM creates a dense and cold shell 
of material consistent with the shell inferred from our results. 

The above results suggest that the high density of the circumburst medium may be 
responsible for the unsuccessful search of an optical counterpart (Jensen et al. 2000, Palazzi 
et al. 2000). In most of the best fit simulations (see Fig. 6), a sizable column density 
of the order of several x 10^^ cm~^ remains after the GRB prompt phase has finished. Dust 
associated to this residual column density may well be responsible for the optical obscuration, 
even for a Galactic dust to gas ratio. In addition, Perna & Lazzati (2002) showed that 
depending on the cloud geometry and GRB spectrum, the dust to gas ratio at the end of 
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lightcurve of the GRB. 
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the burst phase may be even larger than the initial one. 

Remarkably, the slope shown by the X-ray afterglow emission is that typical of a GRB 
with constant density (Sari et al. 1998), which is the assumption for the cloud in the Lazzati 
& Perna (2002) model. We do not find evidence of a temporal break at least until 1 day from 
the main event. Also the afterglow upper limit in T002 is consistent with the same temporal 
slope (see Fig. 5). If, in spite of that, we assume that a break in the light curve occurs 2 
days after the burst, in the scenario of the internal-external shock model, the break could 
be a consequence of a transition to a non relativistic phase of the shock (Dai & Lu 1999), 
due to a high density medium. Using the Dai & Lu (1999) model, we can get an estimate 
of this density. Assuming a Lorentz factor 7 = 1 at t = 2 days, from Dai & Lu (1999), we 
find that E^^ri^ [(1 + z)/2.Qf'/^ = 0.65, where i?54 is the total isotropic-equivalent energy 
E^-ad in units of 10^^ erg, ns is the medium density n in units of 10^ cm~^ and z is the 
GRB redshift. Assuming z = 0.5, which is also consistent with the minimum value inferred 
from the Amati et al. (2002) relation (Ep-Erad) for a 20% dispersion of the data points 
around the best fit, for a standard cosmology with Hq = 71 km s^^ Mpc^^, Qm = 0.3 
and ^A = 0.7 (Spergel et al. 2003), we find E^^ = 0.014 and then a medium density of 
n = 4.5 X 10^ cm^^. This value is somewhat smaller than the best fit densities derived from 
the fit of the column density evolution. However, as shown in the insets of Fig. 6, the data 
do not allow to constrain the cloud geometry with high accuracy, and lower density solutions 
are consistent with the iVn evolution. In addition, given that late afterglow prevents us 
to infer a fading law independently of the late prompt emission time behavior, we cannot 
exclude consistent afterglow solutions with a break before the first NFI observation, thus 
relaxing the upper limit on the ISM density. At any event, for the best fit iVn, the cloud 
geometry and density are consistent with the characteristics of massive cloud cores (Plume 
et al. 1997), and therefore with the progenitor being a massive, short-lived star. 

Continuous, early time monitoring of the time variable opacity (both in single lines and 
in continuum) has the potential to constrain the run of density with radius (Lazzati et al. 
2001b), and hence to trace the mass loss history of the progenitor star. These observations are 
expected to be partly covered by the upcoming satellite Swift, with its repointing capabilities 
of the X-ray and optical detectors in several tens of seconds. A mission which would be 
capable to fully monitor the GRB opacity since the burst onset is the proposed Lobster-ISS 
for the International Space Station, a sensitive all-sky monitor with focusing optics from 
0.1 to 3.5 keV (Fraser et al. 2002) paired with a GRB monitor for higher photon energies 
(from 2 keV to several hundreds of keV) (Frontera et al. 2003b; Amati et al. 2004). When 
measurements of variable absorption in X-rays are complemented with those in the optical 
band, one can further put powerful constraints on the properties of dust in the environment 
of the burst (Perna, Lazzati & Fiore 2003). 
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Table 1. Spectral parameters of the bknpl and CPL best fit models 
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^Tlie column reports the chance probability values obtained with the F-test 
'^The 90% upper boundary value is unconstrained by the fit 
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Table 2. Best fit parameters of the Lazzati & Perna (2002) model for evolution of the 
column density. Error on the parameters are not quoted given the complex covariances (see 

Fig. 6). 

Cloud Geometry Redshift 



sphere 



shell 



dshift 


Nn{0) 


Radius 


x'/dof 




(10^3 cm-2) 


(pc) 




0.1 


6.31 


0.029 


0.65 


0.3 


8.91 


0.086 


0.55 


0.5 


10.0 


0.162 


0.97 


0.1 


3.55 


0.018 


0.15 


0.3 


7.08 


0.036 


0.35 


0.5 


8.91 


0.061 


0.28 
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Fig. 1. — 2-10 keV image of the sky field obtained with MECS in the two BeppoSAX follow- 
up observations of GRB000528. The first TOO started 8.25 hrs after the burst, the second 
3.284 days later. The error box obtained with the WFC 2 is also shown. 
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Fig. 2.— Light curve of GRB000528 in two energy bands, 2-28 keV (WFC), and 40-700 
keV (GRBM), after background subtraction. The zero abscissa corresponds to 2000 May 
28, 08:46:21 UT. The time shces over which the 2-700 keV spectral analysis was performed 
using the 1 s GRBM counters are indicated by vertical dashed hues. The two time shces 
(named 1 and 2 in the text) for which also the 225 channel GRBM spectra were available 
are separated by a vertical arrow. 
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Fig. 3. — Count spectra of the burst in the time intervals A, B, C, D, E, and F. For the 
spectral fits, we summed together the spectra A, B, and C, in order to constrain the model 
parameters. 
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Fig. 4. — Time behaviour of the hydrogen equivalent column density and peak energy of the 
EF{E) spectrum. Also the 2-10 keV and 10-28 keV light curves are shown. 
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Fig. 5. — 2-10 keV fading curve of GRB000528. The afterglow upper limit is at 2a level. 
The vertical dotted line in the Ifet panel shows the beginning of the interval F. Left panel- 
light curve of the GRB and its afterglow assuming as origin the onset time of the GRB. The 
best fit of the light curve in the interval F and its extrapolation at later times is also shown. 
Right panel: The light curve in the interval F and its extrapolation at later epochs, when 
the beginning of the interval F is assumed as origin of the afterglow onset time. 
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Fig. 6. — Redshift corrected column density as a function of the time in the rest frame of the 
burst for three different redshift assumed. Also the best fit model (see text) for the two cloud 
geometries: sphere {continuous line), and shell {dashed line). In the inset the corresponding 
best fit values (initial redshift corrected A^h and cloud radius) are shown for two geometries 
along with the la, 90% and 99% confidence contours. The 90% and the 99% contours are 
not closed. 



